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Abstract 

The demographics of the southern bluefin tuna (SBT) Thunnus maccoyii spawning stock were examined through a large- 
scale monitoring program of the Indonesian longline catch on the spawning ground between 1995 and 2012. The size and 
age structure of the spawning population has undergone significant changes since monitoring began. There has been a 
reduction in the relative abundance of larger/older SBT in the catch since the early 2000s, and a corresponding decrease in 
mean length and age, but there was no evidence of a significant truncation of the age distribution. Pulses of young SBT 
appear in the catches in the early- and mid-2000s and may be the first evidence of increased recruitment into the spawning 
stock since 1995. Fish in these two recruitment pulses were spawned around 1991 and 1997. Size-related variations in sex 
ratio were also observed with female bias for fish less than 170 cm FL and male bias for fish greater than 170 cm FL. This 
trend of increasing proportion of males with size above 170 cm FL is likely to be related to sexual dimorphism in growth 
rates as male length-at-age is greater than that for females after age 10 years. Mean length-at-age of fish aged 8-10 years 
was greater for both males and females on the spawning ground than off the spawning ground, suggesting that size may 
be the dominant factor determining timing of maturation in SBT. In addition to these direct results, the data and samples 
from this program have been central to the assessment and management of this internationally harvested stock. 
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Introduction 

High fishing pressure can affect the size and age structure of a 
fish population, often resulting in a reduction in the relative 
abundance of larger and, indirecdy, older individuals (age-class 
truncation) [1]. In stock assessments, it is often assumed that 
spawning stock biomass (SSB) represents the stock' reproductive 
potential, ignoring changes that may have occurred in the size and 
age structure of the spawning population. Monitoring changes in 
the size and age composition of the spawning population, as well 
as collecting reproductive data to estimate potential annual 
fecundity, is central to quantifying the impacts of fishing on total 
egg production. Direcdy estimating the age structure is particularly 
important for long-lived species where fish size is unlikely be a 
good predictor of age. Changes in the demographics of a spawning 
stock, such as mean length or age and sex ratio, can also be used as 
indicators of the status of the stock for management purposes. 

Southern bluefin tuna (SBT) Thunnus maccoyii is a long-lived, 
late-maturing tuna reaching a maximum size and age of at least 
200 cm and 40 years respectively [2]. SBT has been harvested 
commercially since the early 1950s and is currently estimated to be 
at a small fraction of its original spawning stock size [3]. The 
spawning ground for SBT is located in the north-east Indian 



Ocean between Indonesia and the north-west coast of Australia. 
Mature fish migrate to this region from the Southern Ocean and 
Tasman Sea to spawn between September and April each year 
[4], [5], [6]. An Indonesian longline fishery targeting yellowfin 
[Thunnus albacores) and bigeye tuna (Thunnus obesus), with a bycatch 
of SBT, has operated on the spawning ground since the late- 1970s 
(Table 1) [7], [3]. The catch of SBT by Indonesia increased 
rapidly from just a few tonnes to 2,500 tonnes by 1999 and then 
declined to <650 tonnes in 2003 and 2004. In October 2005, the 
price of fuel in Indonesia doubled after the Indonesian Govern- 
ment cut fuel subsidies. The catch of SBT spiked just prior to the 
price increase as vessels quickly returned to port, but the price rise 
had a significant impact on fishing operations with less longline 
activity and fewer SBT landings [8]. Vessels that continued to 
operate stayed at sea for much longer periods (2-5 months) and 
used fishing vessels as carriers to bring fish back to port [8], [9]. 
Since joining the Commission for the Conservation of Southern 
Bluefin Tuna (CCSBT) in 2008, Indonesia has been allocated a 
national total allowable catch (TAC) ranging between 651 to 750 
tonnes per year. 

Given the significant decline in the SSB and average 
recruitment of SBT [3], it was important to obtain accurate 
estimates of the size and age distribution of SBT landed by the 
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Table 1. Estimated total annual landings (tonnes) of tuna at Benoa, Bali, by species since 1993. 



Year 


Southern bluefin 


Yellowf in 


Bicjeye 


Albacore 


1993 


1 191 


14596 


6192 


716 


1994 


786 


10815 


5360 


338 


1995 


721 


10590 


6543 


463 


1996 


1404 


11061 


10536 


1035 


1997 


1922 


12047 


101 15 


2372 


1998 


1 151 


15735 


1261 1 


905 


1999 


2178 


16128 


9945 


1687 


2000 


1046 


12596 


8676 


2238 


2001 


1419 


12165 


9362 


2461 


2002 


1631 


10380 


1 1646 


2257 


2003 


556 


7399 


5611 


3408 


2004 


613 


4413 


4184 


1906 


2005 


1690 


4196 


3939 


1494 


2006 


558 


4323 


4366 


1450 


2007 


1077 


5354 


5292 


1132 


2008 


905 


6924 


5033 


2811 


2009 


746 


7240 


4680 


1020 


2010 


566 


5372 


2168 


983 


2011 


432 


3006 


2504 


384 


Total 


20592 


1 74340 


1 28763 


29060 



Data sourced from Andaman et al. (2004), Proctor et al. (201 1), and Satria et al. (2012). 
doi:1 0.1 371 /journal.pone.0096392.t001 



Indonesian longline fishery for monitoring changes in the 
spawning population over time and the impact of management 
measures. In 1992, a catch monitoring program was established at 
Benoa (Bali) to examine the catch composition of tunas landed by 
Indonesia [7], [10]. The program focussed on the collection of size 
data of SBT caught by the longline fishery as well as biological 
samples, such as otoliths and ovaries, for analysis of age and 
reproductive dynamics [4], [11], [12]. There were concerns, 
however, that the SBT data from the Indonesian fishery may not 
be representative of the spawning population, given that fish 
caught by Indonesia were generally larger than those caught 
historically by Japan on the spawning ground [13]. Subsequent 
work found that SBT segregate by size on the spawning ground, 
and that this size partitioning was related to spawning activity [1 1]. 
It appears that large fish spend proportionally more time spawning 
while on the spawning ground than small fish, and since spawning 
occurs at the surface, larger fish are more likely to be caught in 
shallow (Indonesian) longline sets, and smaller fish are more likely 
to be caught in deep (Japanese) longline sets. This suggests that, in 
lieu of survey data from which a more complete estimate of the 
size distribution on the spawning ground could be derived, SBT 
caught in shallow longline sets may better represent the spawning 
population. 

This long-term monitoring of the Indonesian fishery and 
collection of large quantities of high quality length and direct 
age data has allowed us to estimate the size and age distribution of 
the Indonesian catch over the past two decades, examine changes 
in size/ age based parameters over time, and estimate sex ratio and 
sex specific growth rates of the spawning population. It has also 
allowed us to examine whether there have been changes in fishing 
practices, such as depth of fishing of the Indonesian fleet, which 



may assist in explaining the changes in the size/age structure of 
the spawning population. 

Methods 

Ethics statement 

Ethical approval was not required for this study, as all fish were 
collected as part of routine fishing procedures. No samples were 
collected by the authors. All samples in this study originated from 
the Indonesian longline fishery and were already dead when 
sampled as part of commercial processing operations. Fish were 
sacrificed by the commercial fisher at sea using standard fisheries 
practices. Permission was granted to use samples from all fish. All 
samples were donated. No field permits were required to collect 
samples, since all originated from commercial catch. SBT are not 
a protected species in any ocean. 

Catch monitoring 

Catch data were obtained from Indonesian longline landings 
through the Benoa-based monitoring program for the 1995 to 
2012 spawning seasons (a spawning season is defined as July 1 of 
the previous year to June 30 of the given year). The data were 
collected through a series of collaborative research programs 
between Australia's CSIRO Marine and Atmospheric Research, 
Indonesia's Research Centre for Capture Fisheries and Research 
Institute for Marine Fisheries, the Indian Ocean Tuna Commis- 
sion, and Japan's Overseas Fishery Cooperation Foundation [8]. 
In mid-2002, the monitoring program expanded to include the 
ports of Muara Baru (Jakarta) and Cilacap (south coast Central 
Java); however, the majority of targeted SBT sampling still 
occured at Benoa as this is the port where the majority (85%) of 
SBT is landed [14]. Landed SBT are graded into export and non- 
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2007 



2008 



2009 



2010 



2011 



20 25 30 35 40 
Age (years) 



170 190 
Fork length {cm 



Figure 1 . Length (left) and age (right) distribution of the SBT spawning stock by spawning season. The dark grey bar shows the median 
length class. For comparison, the length distribution of SBT thought to be caught south of the spawning ground (Company A) is shown for the 2004 
(n = 121), 2005 (n = 685), 2006 (n = 31 1) and 2007 (n = 452) seasons (solid grey line). The 3-year running mean of the age distribution is shown (solid 
black line) and the two possible pulses of recruitment to the fishery (dashed grey line). 
doi:1 0.1 371 /journal.pone.0096392.g001 



export quality based on flesh quality. Quality and grading is 
dependent on handling, length of trip and/ or condition of fish at 
capture, rather than fish size [1 1]. Generally, only SBT graded as 
not suitable for export were available for monitoring. 



Fork length (FL) was measured to the nearest cm (Table 2) and 
length frequencies were constructed for each spawning season. 
Through the monitoring program, one fishing company was 
identified as having vessels operating well south of the spawning 
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Figure 2. Proportion of SBT caught by size (A) and age (B) classes by spawning season. Mean fork length (FL) ± 2SE (upper graph) and 
mean age of all and SBT >20 years (yrs) (lower graph) by spawning seasons are shown in red. Data from Company A operating south of the spawning 
ground were excluded. 
doi:1 0.1 371 /journal.pone.0096392.g002 



ground (southern zone) from January 2004 to April 2007. Since it 
is unknown if fish caught in the southern zone were mature or 
would migrate to the spawning ground, it was important that these 
fish were identified in the size data and excluded from our analysis 
so that the estimated size distribution of the spawning population 
compared over time was consistent. The size distribution of SBT 
landed by this company (Company A) was compared to the other 
companies (combined data) each season using the two-sample 
Kolmogorov-Smirnov (K-S) test to determine if significant 
differences were present. 

Size distribution 

To examine changes in the size structure of the spawning 
population over time, it was necessary to determine whether 
changes had occurred in the fishing practices of the Indonesian 
fleet, such as the depth of the longline sets or distribution of effort 
through the spawning season, which may influence proportion of 
the population available to be caught and the size offish caught. A 
continuous, consistent time series of data on fishing depth are not 



available for the Indonesian longline fishery. As an alternative we 
used the catch data obtained through the monitoring program to 
calculate the relative proportion of bigeye tuna in the landings, 
referred to as the bigeye (BE) index [11]: BE index = weight of 
bigeye/ (weight of bigeye + yellowfin). 

The index is based on the assumption/ observation that bigeye is 
generally caught deeper than yellowfin. The BE index was 
calculated for individual landings from 1995-2009, and was 
available for 96.8% of the landings in those years which had SBT 
length measurements. The BE index ranges from 0 to 1 and is 
divided into 5 levels (0.0-0.2, 0.2-0.4, etc.) for analysis. If the BE 
index is high, the fishing depth is assumed to be deep (and vice 
versa). To take into account possible differences in the size 
distribution of SBT at different depth levels (i.e., different levels of 
the BE index), we calculated the weighted mean length of SBT 
caught each spawning season using weights equal to the 
proportion of fish at each BE index level. Similarly, to examine 
the effect of month of capture on the size distribution of the catch, 
we calculated the weighted mean length of SBT caught each 
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1994 1996 1998 



2000 2002 2004 
Spawning season 



2006 2008 2010 



Figure 3. Mean bigeye index (x) and proportion <165 cm FL (•) by spawning season. 

doi:1 0.1 371 /journal.pone.0096392.g003 



spawning season using weights equal to the proportion of fish 
caught in each month. The only months that had length data 
collected in all seasons were September to March, thus length data 
for April to August were excluded (0.3% of length measurements). 
A comparison was then made with the mean fish length by 
spawning season calculated using the original sample data (i.e., the 
unweighted means). 

To formalize the above investigations, generalized linear models 
(GLMs) were fit to the data, with the proportion of small SBT (< 
165 cm FL) as the binomial response and spawning season, BE 
index and month as potential explanatory variables. The variables 
and interactions terms to include in the final model were 
determined using % 2 significance tests. Using the best fitting 
model, the expected proportion of small SBT in each season was 
estimated for fixed levels of month and BE index. 

Age estimation 

Sagittal otoliths were removed from approximately 500 to 1 700 
SBT measured for length each spawning season (1995-201 1), with 
the exception of the 1996 season when only 225 otoliths were 
sampled (Table 2). Five hundred were selected from each season 
for annual age estimation apart from 1996 when only 50 were 
selected. The number selected was based on the work by Morton 
and Bravington [15] who estimated that 500 samples would be 
sufficient to provide acceptable precision; i.e., coefficients of 
variation (CVs) under 20% for the Indonesian fishery. For the 
1995 to 1999 seasons, otoliths were selected randomly from those 
sampled (which included most otoliths collected). For the 
remaining spawning seasons, a fixed number of otoliths were 
chosen from each 1 cm length class to obtain as many age 
estimates from length classes where sample sizes were small. 

All otoliths were prepared and read following the techniques 
described by Clear et aJL [16] and Anonymous [17] and the 
precision was examined by calculating the CV of replicates 
readings [18]. To determine the age structure of the Indonesian 
catch, age-length keys (ALK) were developed using the sample of 
aged fish for each spawning season except for the 1996 season as 
too few otoliths were read to develop a reliable ALK. The ALK 
gives the proportion offish at age in each 5-cm length class, which 
is then used to estimate the age-frequency distribution of the catch 
from the length-frequency distribution obtained through the 



monitoring. The Shannon-Weaver index was calculated as a 
measure of age structure diversity, using: 

H' = -J2pM.(Pi) 

where p, is the proportion offish belonging to the z'th age class [19]. 

Sex ratio 

Sex was recorded for the majority offish measured for length for 
the 2000-2012 spawning seasons. For most fish, sex was 
determined from a small sample of remnant gonad (RG) tissue 
left in the visceral cavity when the fish was landed, as SBT lack 
external sexual characteristics and almost all were landed in gutted 
form. A secondary method, based on the size and shape of the 
anus, was used to identify sex when RG tissue was not present. 
The method used to determine sex was recorded only for the 
2009-2012 seasons, thus it was unknown whether the RG or 
'anus' method was used in the earlier seasons. Exploratory data 
analysis suggested that the 'anus' method was not accurately 
identifying the sex of all fish, thus only RG data for 2009-2012 
were used here (Table 2). The sex ratio of SBT was calculated, and 
chi square tests were used to examine differences from an expected 
1:1 by season and length class (5-cm). 

Sex-specific growth rates 

The distribution of size-at-age between males and females on 
the spawning grounds were examined using samples for which sex 
was determined by the RG method (from seasons 2009-2012). A 
von Bertalanffy (VB) growth curve was fit to the age and length 
data for each sex separately. Although SBT growth appears to 
follow a two-stanza VB model, the transition between growth 
stanzas occurs before age 5 [20]. Since the age of transition is 
younger than any data from the spawning ground, fitting a simple 
VB model was considered the most parsimonious approach. 

In addition, age and length data for SBT caught off the 
spawning ground (see [12]) were used to compare length-at-age of 
SBT on and off the spawning ground. Only samples with sex 
information were used so that males and females could be 
compared separately. For all fish caught off the spawning ground, 
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0.6 




1993 1995 1997 1999 2001 2003 2005 2007 2009 2011 2013 

Spawning season 

Figure 4. Proportion of SBT <1 65 cm by spawning season for bigeye index (A) and month-group (B). BEI0 is shallow, BEI4 is deep. BEI = 
bigeye index. 

doi:10.1371/journal.pone.0096392.g004 



sex was determined by trained observers onboard vessels based on 
the appearance of the whole gonad [12]. 

Results 

Length distribution 

Length measurements were obtained for nearly 23,000 SBT 
from 18 spawning seasons (Table 2). Length data from the landed 
catch of Company A, identified as having fishing vessels operating 
in the southern zone, was significandy different to the other 
companies in the 2005 to 2007 seasons (K-S test; P<0.001). The 
catch by Company A clearly comprised a greater proportion of 
small fish (140-160 cm FL) in each of the seasons in question 
(Fig. 1). The catch by Company A in the 2004 season also 



appeared to contain more small fish relative to the other 
companies, but the difference was not significant (K-S test; 
P— 0.053). Since it was unknown if small SBT caught in the 
southern zone would migrate to the spawning ground, SBT landed 
by Company A in the 2004 to 2007 seasons were excluded from 
further analyses. 

Considerable change has occurred in the size distribution of 
SBT caught by Indonesian longliners on the spawning ground 
since monitoring began. In the mid- to late- 1990s, the majority of 
fish caught were between 165 and 190 cm FL with a median 
length of approximately 180 cm FL (Fig. 1). In 1998, an increase 
in the relative proportion of small fish (< 1 65 cm FL) occurred but 
this did not persist into the 1999 season (Fig. 2A). In 2000, the 
relative proportion of small fish increased again compared to the 
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184 




166 H 1 1 1 1 — 

1992 1997 2002 2007 2012 



Spawning season 

Figure 5. Mean fork length of SBT by spawning season. Data are shown for observed (A) and estimated mean length when BE index (•) and 
month (x) is accounted for in the length distribution of fish caught. 
doi:1 0.1 371 /journal. pone.0096392.g005 

previous season, and continued to increase steadily, peaking at 
45.9% of the catch in 2003 (Fig. 2A). Since then, the relative 
abundance offish < 1 65 cm FL has fluctuated between 22.5 and 
36.9% with a median length of approximately 170 cm FL. The 
mean length of SBT in the catch decreased between 1995 and 
2003 from 180.7 cm to 166.8 cm FL, and has remained between 
168.3 and 171.0 cm FL since that time (Table 2; Fig. 2A). 

Initial examination of mean BE index by spawning season 
suggests that there has been no systematic change in fishing depth 
by the Indonesian fishery that would account for the increase in 
the relative abundance of small SBT in the 2000s (Table 2). 
Although there appears to be a relationship between mean BE 
index (fishing depth) and the proportion of fish < 1 65 cm FL in the 
catch between 1995 and 1999 (Fig. 3), this relationship was not 
present in the latter years when the greatest change in the size of 
fish caught occurred. The increase in the proportion of fish < 
165 cm FL in the catch from 2000 occurred at all levels of the BE 
index and month groups (Fig. 4) suggesting that it occurred 
independently of changes in fishing practices. This is supported by 
the comparison of estimated mean length by spawning season 
when BE index and month of capture were accounted for (Fig. 5). 
The similarity of the observed and estimated weighted mean 
length of fish by season suggests that the observed changes in 
fishing practices had only a minor influence on the size of fish 
caught. 

Results from the GLMs suggest that spawning season, month 
and BE index are all highly significant explanatory variables of the 
proportion of small SBT, as are the pair-wise interaction terms 
between spawning season and month and spawning season and BE 
index (Table 3). The interaction term between month and BE 
index did not significantly improve the model fit (x 2 test, 
P— 0.235), so was not included in the final model. Using the final 
model to estimate the expected proportion of small SBT at each 
level of month and BE index confirms that the proportion of small 
SBT on the spawning ground increased between seasons 2000 and 



2003 (Fig. 6). Although the size of this increase varied among some 
months and levels of the BE index it was consistently present. 

Age distribution 

Age was estimated for 7,773 of the 8,000 SBT selected for age 
estimation. These fish ranged in size from 122 to 229 cm FL and 5 
to 40 years. The CV between readings by the primary reader 
ranged from 3.5-5.3% by spawning season with an overall CV of 
4.15%. The second age estimate of the primary reader agreed with 
the original estimate in 37.6% of cases, and was within two years 
of the original in 91.7% of cases. The CV between primary and 
secondary readers was 5.08%. These low levels of error, especially 
between the two readers, suggest consistent interpretation of age in 
blind tests. 

The minimum age of SBT sampled was 5 years for two unsexed 
fish, while the maximum age for males and females was 38 years 
and 40 years, respectively. As expected, the estimated age 
structure of fish caught by the Indonesian fleet on the spawning 
ground shows substantial changes over time (Fig. 1; Fig. 2B). 
Between 1995 and 2000, the majority offish landed were relatively 
old (&20 years) and the mean age was 19-21 years (Fig. 2B; 
Table 2). After 2000, the mean age offish decreased to 14-17 
years as the abundance of young fish increased in the catch relative 
to old fish (Fig. 1; Table 2). The age structure diversity index was 
significantly higher (ANOVA, F= 184.5, P<0.001) in the 1998- 
2000 period (mean 22.08±1.3) compared to the 2003-2011 
period (mean 16.3±0.55). A clear increase in the relative 
abundance of young fish occurred in the early 2000s and again 
in ~2006, and these modes appeared to progress through the 
fishery on an annual time step (Fig. 1). The first pulse has a clear 
mode in 2003 (age 12 years) which follows through to 201 1 (age 20 
years). These fish would have been spawned in 199 1 . Similarly, the 
second pulse has a clear mode in 2006 (age 9 years) which follows 
through to 2011 (age 14 years). These fish would have been 
spawned in 1997. The average age of SBT greater than 20 years 
old has remained relatively stable since monitoring began, 
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Figure 6. Expected proportion of small SBT each month and level of the bigeye index. Proportions calculated from the best-fitting 
generalized linear model. BEIO is shallow, BEI4 is deep. BEI = bigeye index. 
doi:1 0.1 371 /journal.pone.0096392.g006 



although there is some indication of a slight decline in the most 
recent four or five spawning seasons (Fig. 2B). 

Sex ratio 

The sex ratio of SBT in the landed catches was significantly 
different from the expected 1:1 in 2009 to 2011 with a slight 
dominance of females (x 2 tests, P<0.01) (Table 2). Although the 
percent females was slightly higher in 2012 than in 2009 to 201 1, 
the sex ratio was not significantly different from 1:1 in 20 1 2 due to 
the smaller sample size (x 2 test, P= 0.066). The sex ratio varied 
with fish length in all seasons, being significantly female-biased in 
length classes below 170 cm FL, and significantly male-biased in 
length classes above 170 cm FL (x 2 tests, P<0.01) (Fig. 7). Chi- 
square tests did not indicate significant deviations from a 1 : 1 ratio 
in the 1 70 cm length class in any season. A clear trend was evident 
in the percent female with fish length, decreasing steadily from 74— 
89% in the 150 and 155 cm length classes to only 7-27% in the 
185 and 190 cm length classes. 

Sex specific growth rates 

Mean length-at-age for males on the spawning grounds is 
significantly greater than that of females beyond age 10 (Fig. 8). 
VB growth models fit to the data suggest males grow at a faster 
rate than females, but that females have a higher asymptotic length 
and, thus, 'catch up' in size by age 35 (Fig. 9). The higher 
asymptotic length for females may be due to small sample sizes 
beyond age 25 (see Discussion). Comparing fish caught on and off 
the spawning grounds show that the mean length-at-age of fish 
aged 8-10 years is significantly greater on the spawning ground, 
for both males and females (Fig. 10). 

Discussion 

Our results clearly show that the size and age distribution of 
SBT caught in the Indonesian tropical tuna longline fishery on the 
SBT spawning ground has undergone substantial changes over the 
past two decades. These changes appear to be irrespective of 
changes in fishing practices for which we have information, such as 
fishing depth (BE index) or month of capture. The catch 
monitoring showed that there has been a decrease in the relative 
abundance of larger/older SBT since the early-2000s. It was 
important, however, to identify fish in the size monitoring data 



that were caught south of the spawning ground so that these could 
be excluded from the analyses. SBT catches landed by the 
company identified as having operated in this south zone 
contained a greater proportion of small fish (140-160 cm FL) 
compared to the other companies in the 2005 to 2007 seasons and 
to a lesser extent in 2004. This is consistent with historic Japanese 
catch data which showed that the mode of SBT caught on the 
spawning ('Oka') ground at 10-20°S was higher than for the 
staging ('Old') ground to the south at 20-35°S [21]. Although the 
maturity status of small fish caught in the southern zone was 
unknown, it is possible that a proportion were immature and 
undertaking trial migrations towards the spawning ground. 
Biological sampling and histological analysis of gonads is required 
to confirm the maturity status of the smaller 'southern' fish. 

Changes in the size of SBT caught on the spawning ground 
have been reported previously. The average size of SBT caught by 
Japanese longliners increased steadily between the early- 1970s 
(161.0 cm FL) and early-1990s (169.9 cm FL), and it was suggested 
that the increase was due to a reduction in the number of small 
SBT in the spawning population or an increase in growth rates 
over time [13]. In the late- 1990s, the mean size of SBT caught on 
the spawning ground during longline surveys by Japan was lower 
than that caught by Indonesia, and these differences were 
explained by different fishing depths targeted by the vessels [11]. 
The decrease in the mean size (and age) of SBT caught by 
Indonesia since the early 2000s appears to be independent of 
changes in fishing depth, but should be considered in relation to 
effort within the fishery to determine if the change was the result of 
an increase in number of small/young fish or a decrease in the 
larger/old fish in the catch. Up until 1997, it was possible to 
determine catch per unit effort (CPUE) for most Indonesian vessels 
because it was known how many days vessels had fished and the 
number of hooks used because the crew were interviewed when 
they landed their catch. Carrier boats were first noted in the 
fishery in 1997 when four were registered as using Benoa [22]. 
This information on the number of sets and the number of hooks 
used per set was lost as carrier boats became increasingly 
dominant in the fishery and crews of individual boats could not 
be questioned as they did not land their catch. Despite the lack of 
CPUE information, the decrease in mean length and age of SBT 
observed in the current monitoring data appears to be due, in part, 
to pulses of higher than average abundance of young fish entering 



Table 3. Analysis of deviance table from fitting a generalized linear model to the proportion of small (<165 cm) SBT on the 
spawning ground. 







Df 


Deviance 


Resid. Df 


Resld. Dev 


PO|Chi|) 


NULL 






15814 


21923 




Season 


14 


6388.8 


15800 


15534 


0.000 


Month 


6 


154.9 
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15379 


0.000 


BE Index 


4 


37.6 


15790 


15342 


0.000 


Season:Month 


78 


302.9 


15712 


15039 


0.000 


Season:BE Index 


52 


133.8 


15660 


14905 


0.000 


Month:BE Index 


24 


28.6 


15636 


14876 


0.235 



Results are for the model including spawning season, month and bigeye index as main effects and all pair-wise interactions. 
doi:1 0.1 371 /journal.pone.0096392.t003 
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Figure 7. Percent female by 5-cm length class for SBT in the Indonesian catch in 2009-2012. Data point excluded if n- 10. 
doi:1 0.1 371 /journal.pone.0096392.g007 



the fishery in the early- and mid-2000s. The absence of a decline 
in the mean age offish &20 years in the early 2000s suggests that 
the shift in the age distribution of SBT caught was unlikely to be 
due to a decrease in the catch of very old fish. The slight decrease 
in mean age of fish &20 years in 2008-201 1 may be attributed to 
the first pulse of recruits entering the 20+ age class in these recent 
years. The two clear pulses of recruitment in the age frequency 



data were not apparent in the length frequency time series. This 
can be explained by the substantial variability in length-at-age of 
adult SBT obscuring these cohorts in the length frequency data 
[2], [12]. _ 

It is possible that the observed increase in the relative number of 
smaller/younger fish in catch since the early 2000s occurred as a 
result of increased juvenile survival following the introduction of 
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quotas in the mid- 1 980s, leading to more fish surviving though to 
the spawning age. Quotas have reduced the total catch of SBT 
from around 40,000 tonnes in the early- 1980s to around 10,000 
tonnes since 2006 [23]. Importantly, the introduction of quotas 
significandy reduced the number of juveniles (aged 2-4 years) 
being caught in surface fisheries around Australia from a peak of 
just over 21,000 tonnes in 1982 to around 5,000 tonnes since 
1990, allowing them to escape to the high-seas longline fisheries. 
However, the pulse offish entering the spawning population as say 
12 years-olds in 2003 would have been 2-4 year-olds in 1993- 
1995, after quotas had taken effect, suggesting that both natural 
variability in recruitment strength and the effect of fishing shape 
the age structure of the SBT spawning population. Recruitment of 
young SBT to the New Zealand longline charter fishery is evident 
in the early-2000s and more recendy in 2008 [24]. The later pulse, 
aged 6-8 years in 2012, may become evident as the beginning of a 
new pulse of increased abundance in the spawning stock at age 9 
years in 2014. 

The catch of SBT by Indonesia' longline fishery was dominated 
by females in the smallest length classes (145-165 cm FL), while 
the catch of SBT by Japanese longliners on feeding grounds in the 
southern oceans were close to 1:1 for these length classes [12]. It 
could be suggested that the determination of sex may not be 
accurate for SBT in Indonesia since it was done using a small 
sample of remnant gonad tissue left in the visceral cavity when the 
fish is landed. Although this method is crude and there is potential 
for error, the length-at-age data indicated clear sexual dimorphism 
in growth suggesting that sex has been correctiy identified in most 
cases. The shift to male dominance in the current study occurred 
at approximately the same length as observed in [12], and the 
trend of increasing male dominance with size was also very similar. 
Sex ratios reported for other tuna species are generally 1 : 1 with a 



similar prevalence of males in the larger length classes [25]. The 
male dominance in large length classes for SBT may be due to the 
observed sexual dimorphism in growth, as suggested for albacore 
tuna (Thumus alalungd) in the South Pacific Ocean [26]. 

The reason for the predominance of females in length classes < 
170 cm is less clear but may reflect sexual differences in 
vulnerability or availability on the spawning ground. Females 
may be more catchable than males if they are more likely to be 
feeding or actively spawning (i.e., spending a greater proportion of 
time near the surface targeted by the shallow-setting Indonesian 
longliners) compared to males of a similar size. Unfortunately, 
data on feeding, depth partitioning, residency and spawning 
behaviour are not available for SBT by sex, although females are 
known to be capable of spawning daily [4]. Regardless of the cause 
of the bias in sex ratio, a disproportional harvest of females is likely 
to have implications for population egg production estimates, 
modeling stock dynamics and fisheries management. A higher 
catch rate of females by Indonesia over time could ultimately lead 
to a decrease in the abundance of females in the spawning 
population, and a subsequent decline in reproductive potential. 

The growth rate of male SBT on the spawning ground was 
found to be greater than that of females, supporting results in Lin 
and Tzeng [27]. However, unlike Lin and Tzeng [27], our results 
show a higher asymptotic length for females. This could be due to 
small sample sizes at the oldest age, because we did not find it to be 
true when we refit VB models to all of the data, regardless of 
sexing method (in which case asymptotic length was estimated to 
be equal for males and females). Sexual dimorphism in the growth 
of other tuna species has also been found (e.g. albacore [28]). 
Mean length-at-age of fish aged 8—10 years, both males and 
females, was greater on the spawning ground than off. This is 
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Figure 10. Mean length-at-age ± 2SE for fish caught on and off the spawning ground by sex. SG = spawning ground. 
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consistent with [29] which showed that size is a factor in the early 
maturation of SBT for both sexes. 

The observed changes in the length and age structure of the 
catches since the mid-1990s may not reflect a reduction in the per 
capita reproductive potential of the spawning stock. Although 
there has been a decrease in the mean size and age of adults 
caught by Indonesia on the spawning ground, the change in 



demographics may be due to an increased abundance of young 
fish recruiting to the fishery rather than a substantial truncation of 
the length and age distribution due to fishing mortality. The first 
pulse of recruits that appeared in the catch in the early 2000s are 
now reaching ages of ~20 years and their higher batch fecundity 
relative to younger females [4] may mean that per capita annual 
egg production could increase over time. 
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The size- and age-based metrics presented here are essential to 
understanding the impact of the Indonesian fishery on the SBT 
spawning population and form part of a suite of indicators on the 
status of the SBT stock [30]. The age estimates are used to 
estimate catch-at-age for the Indonesian fishery which is an 
important input to the stock assessment models providing 
information on the adult component of the population and 
changes over time [31]. These data are also used in the CCSBT 
operating model which has been used for Management Strategy 
Evaluation of the adopted SBT management procedure. The 
management procedure is used to set the global TAC of SBT [32]. 
The continued collection of the time series of data will be vital for 
monitoring the rebuilding of the stock and the international 
management of the species. The monitoring program and these 
data are also an important component of a new ground breaking 
method for estimating the abundance of the spawning population 
of SBT using modern genetic techniques to identify parent- 
offspring-pairs - the "close-kin" method [33]. However, uncer- 
tainties remain in estimating the reproductive potential of SBT. 
Data is required to accurately estimate the maturity schedule for 
SBT, understand spawning migrations and reproductive behavior, 
as well as determine residency and selectivity on the spawning 
ground. 

References 

1 . Berkeley SA, Hixon MA, Larson RJ, Love MS (2004) Fisheries sustainability via 
protection of age structure and spatial distribution of fish populations. Fisheries 
29: 23-32. 

2. Gunn JS, Clear NP, Carter TI, Rccs AJ, Stanley CA, ct al. (2008) Age and 
growth in southern bluefin tuna, Thunnus maccoyii (Castelnau): direct estimation 
from otoliths, scales and vertebrae. Fish Res 92: 207-220. 

3. CCSBT (2013a) Report of the eighteenth meeting of the Scientific Committee. 
Commission for the Conservation of Southern Bluefin, 7 September, Canberra, 
Australia. Available: http://www.ccsbt.org/site/reports_past_meetings.php. Ac- 
cessed 2013 Dec 2. 

4. Farley JH, Davis FLO (1998) Reproductive dynamics of southern bluefm tuna, 
Thunnus maccoyii. Fish Bull 96: 223-236. 

5. Evans K, Patterson TA, Reid H, Harley SJ (2012) Reproductive schedules in 
southern bluefin tuna: Are current assumptions appropriate? PLoS ONE 7(4): 
c34550. doi:10.1371/journal.pone.0034550. 

6. Patterson TP, Evans K, Carter TI, Gunn JS (2008) Movement and behaviour of 
large southern bluefin tuna {Thunnus rnaccciyii) in the Australian region determined 
using pop-up satellite archival tags. Fish Oceanog 17: 352—367. 

7. Davis TLO, Bahar S, Farley JH (1995) Southern bluefin tuna in the Indonesian 
longiine fishery: historical development, composition, season, some biological 
parameters, landing estimation and catch statistics for 1993. Indones Fish Res J 
1: 68-86. 

8. Proctor CH, Andamari R, Retnowati D, Herrera M, Poisson F, et al. (2006) The 
catch of SBT by the Indonesian longiine fishery operating out of Benoa, Bali in 
2005. Report CCSBT-ESC/0609/10 prepared for the 11 th CCSBT Extended 
Scientific Meeting 4-1 1 September, Tokyo, Japan. 

9. Proctor C, Purwanto, Budhiman A, Nurhakim S, Prisantoso BI, et al. (2011) 
Project Summary - ACIAR Project FIS/2002/074 Capacity development to 
monitor, analyse and report on Indonesian tuna fisheries. Australian Centre for 
International Agricultural Research, Canberra. 19 p. Available: http://aciar. 
gov.au/publication/cop017. Accessed 2013 Dec 2. 

10. Davis TLO, Bahar S, Naamin N, Farley JH, Lc D (1999) Trends in tuna and 
billlish landings by the longiine fishery operating out of Benoa, Bali. Indones 
Fish Res J 5: 50-60. 

11. Davis TLO, Farley JH (2001) Size partitioning by depth of southern bluefin 
tuna, Thunnus maccoyii, on the spawning grounds. Fish Bull 99: 381—386. 

12. Farley JH, Davis TLO, Gunn JS, Clear NP, Preecc AL (2007) Demographic 
patterns of southern bluefin tuna, Thunnus maccoyii, as inferred from direct age 
data. Fish Res 83: 151-161. doi:10.1016/j.fishres.2006.09.006. 

13. Gunn JS, Davis TLO, Farley JH, Clear NP, Haskard K (1996) Preliminary 
estimations of the age structure of the SBT spawning stock. Report CCSBT/ 
SC/96/10 prepared for the 2 nd CCSBT Scientific Meeting 26 Aug - 5 Sep, 
Hobart, Australia. 

14. Satria F, Suman A, Nugroho D, Nugraha B, Widodo A, ct al. (2012) Indonesia 
southern bluefin tuna fisheries, a national report year 2011. Report CCSBT- 
ESC/1208/SBT Fisheries - Indonesia prepared for the 17 th CCSBT Meeting of 
the Scientific Meeting 27—31 August, Tokyo, Japan. 

15. Morton R, Bravington M (2003) Estimation of age profiles of southern bluefin 
tuna. Report CCSBT-ESC/0309/32 prepared for the CCSBT Scientific 
Meeting 1-4 September, Christchurch, New Zealand. 



Acknowledgments 

The success of the SBT monitoring program was only possible due to the 
dedicated efforts of all the scientist team at Research Institute for Tuna 
Fisheries (Benoa, Bali), and in particular that of K. Siregar and R. Mashar, 
the biological samplers and primary measurers of the SBT. We also thank 
E. Rahmat (RIMF) for earlier years of data entry into the SBT biological 
database and to B. Iskandar Prisantoso (RCFMC) in his former role as 
database manager. The cooperation of the longiine tuna industry- 
coordinated through Asosiasi Tuna Longiine Indonesia) and the individual 
processing companies in providing access and facilities to carry out the 
sampling is much appreciated. We also acknowledge the support of all 
other agencies within Ministry of Marine Affairs and Fisheries (Indonesia) 
for the research activities. In addition, wc thank T. Carter for cleaning and 
archiving the otoliths. This manuscript was improved by comments from 
A. Preece and the anonymous reviewers. 

Author Contributions 

Conceived and designed the experiments: JHF TLOD CHP CRD. 
Analyzed the data: JHF JPE TLOD. Contributed reagents/materials/ 
analysis tools: JHF JPE TLOD RA CHP BN CRD. Wrote the paper: JHF 
JPE TLOD RA CHP BN CRD. 



16. Clear NP, Gunn JS, Rees AJ (2000) Direct validation of annual increments in the 
otoliths of juvenile southern bluefin tuna, Thunnus maccoyii, through a large-scale 
mark-and-rccapturc experiment using strontium chloride, lush Bull 98: 25-40. 

17. Anonymous (2002) A manual for age determination of southern bluefin tuna 
Thunnus maccoyii - otolith sampling, preparation and interpretation. The direct 
age estimation workshop of the CCSBT, 11-14 June, 2002, Queenscliff, 
Australia. Available: http:/ /www. ccsbt.org/userf lies/ file/ docs_english/ 
mectings/meeting_reports/ccsbt_09/report_of_daews.pdf. Accessed 2013 Dec 
2. 

18. Campana SE, Anmand MC, McMillan JI (1995) Graphical and statistical 
methods for determining the consistency of age determinations. Trans Am Fish 
Soc 124: 131-138. 

19. Shannon CE (1948) A mathematical theory of communication. Bell Syst TechJ 
27: 379-423, 623-656. 

20. Polacheck T, EvesonJP, Laslett GM (2004) Increase in growth rates of southern 
bluefin tuna {Thunnus maccoyii) over four decades: 1960 to 2000. CanJ Fish Aquat 
Sci 61(2): 307-322. 

21. Shingu C (1978) Ecology and stock of southern bluefin tuna. Fish Study Jap 
Assoc Fish Resources Protection, 31: 81 p. (In Japanese; English translation in 
Rep CSIRO Div Fish Oceanogr 131: 79 p. 1981.) 

22. Davis TLO, Bahar S, Naamin N, Lc D (1998) Catch monitoring of the fresh 
tuna caught by the Bali-based longiine fishery. Report CCSBT/SC/98/5 
prepared for the 4 th CCSBT Scientific Meeting 23-31 July, Shimizu, Japan. 

23. CCSBT (2013b) Estimated Total Global Catch of Southern Bluefin Tuna. 
Available: http://www.ccsbt.org/site/sbt_data.php. Accessed 2013 Dec 2. 

24. Anonymous (2013) Annual review of national SBT fisheries for the Scientific 
Committee. Report CCSBT-ESC/1309/SBT Fisheries - New Zealand prepared 
for the 18 ,h Meeting of the CCSBT Scientific Committee, 2-7 January, 
Canberra, Australia. 30 p. 

25. Schaefer KM (2001) Reproductive biology of tunas. In: Block BA, Stevens ED, 
editors. Tuna: physiology, ecology and evolution, Vol 19. San Diego: Academic 
Prcss.pp. 225-270. 

26. Farley JH, Williams AJ, Hoyle SD, Davies CR, Nicol SJ (2013) Reproductive 
dynamics and potential annual fecundity of South Pacific albacore tuna {Thunnus 
akilunga). PLoS ONE 8(4): c60577. doi:10.1371/journal.pone.0060577. 

27. Lin Y-T, Tzeng W-N (2010) Sexual dimorphism in the growth rate of southern 
bluefin tuna Thunnus maccoyii in the Indian Ocean. J Fish Soc Taiwan 37(2): 135- 
151. 

28. Williams AJ, Farley JH, Hoyle SD, Davies CR, Nicol SJ (2012) Spatial and sex- 
specific variation in growth of albacore tuna {'Thunnus alalunga) across the South 
Pacific Ocean. PLoS ONE 7(6): c.39318. doi:10.1371/journal.pone.0039318. 

29. Farley JH, Davis TLO, EvesonJP (2001) Length and age distribution of SBT in 
the Indonesian longiine catch on the spawning ground. Report CCSBT/SC/ 
0108/12 prepared fbr the 2nd CCSBT Stock Assessment Meeting 19-28 
August, Tokyo, Japan. 

30. Patterson H, Preece A, Hartog J (2012) Fishery indicators for the southern 
bluefin tuna stock 2011-12. Report CCSBT-ESC/1208/14 prepared for the 
17 th Meeting of the CCSBT Scientific Committee, 27-31 August, Tokyo, Japan. 

31. Hillary RM, Preece A, Davies C, Barnes B, Bcgg G (201 1) Reconditioning of the 
southern bluefin tuna operating model: exploratory data analysis, fitting 



PLOS ONE | www.plosone.org 



14 



May 2014 | Volume 9 | Issue 5 | e96392 



Southern Bluefin Tuna Spawning Ground Demographics 



performance, and current stock status. Report CCSBT-ESG/1 107/1 1 prepared 
for the 16 th CCSBT Meeting of the Scientific Meeting 19-28 July, Bali, 
Indonesia. 

32. Hillary R, Preece A, Davies C (2012) Developing a management procedure 
based recovery plan for southern bluefin tuna. Final report for FRDG project 



2011/034. Hobart: CSIRO. 131 p. Available: http://frdc.com.au/research/ 
final-reports/Pages/defaultaspx. Accessed 2014 Apr 8. 
33. Bravington MV, Grewe PG, Davies CR (2014) Fishery-independent estimate of 
spawning biomass of southern bluefin tuna through identification of elose-kin 
using genetic markers. Final report for FRDC project 2007/034. Hobart: 
CSIRO. 147 p. 



PLOS ONE | www.plosone.org 



15 



May 2014 | Volume 9 | Issue 5 | e96392 



